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ABSTRACT 

We measured the oxygen abundances of the warm (T~ 10 4 K) phase of gas in 
seven early-type galaxies through long-slit observations. A template spectra was 
constructed from galaxies void of warm gas and subtracted from the emission- 
line galaxies, allowing for a clean measurement of the nebular lines. The ratios 
of the emission lines are consistent with photoionization, which likely originates 
from the UV flux of post-asymototic giant branch (PAGB) stars. We employ 
H II region photoionization models to determine a mean oxygen metallicity of 
1.01 ± 0.50 solar for the warm interstellar medium (ISM) in this sample. This 
warm ISM 0.5 to 1.5 solar metallicity is consistent with modern determinations of 
the metallicity in the hot (T~ 10 6 - 10 7 K) ISM and the upper range of this warm 
ISM metallicity is consistent with stellar population metallicity determinations. 
A solar metallicity of the warm ISM favors an internal origin for the warm ISM 
such as AGB mass loss within the galaxy. 



Subject headings: Galaxies: Elliptical and Lenticular, Galaxies:ISM, Galaxies:Abundances 
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Introduction 



1.1. Warm ISM Discovery 



Early-type gala xies were once thought t o contain 



very 



li ttle gas in their interstell ar 



medium (ISM) (e.g. Mathews k Baked (119711 ) : lBregmanl (119781 ): IWhite fc Chevalied (119831 ^. 
but are now kn own to host massive amounts of extremely hot (T~ 10 6 — lO 7 -^) gas that emits 



in the X-rays (jForman et al.lll979l ). The reason early- type galaxies were originally thought 
to be void of gas is that photographic-plat e surveys did not reveal the warm (T~ lO 4 ^ ) 
ISM that is abundant in spiral galaxies (e.g. iMayalll (119361 ); ISandagd (119571 ); iMayalll (119581 )). 
Only with the advent of large telescopes and efficient detectors was the warm ISM detected, 
although the masses are an order of magnitude lower than that of spirals. Spectroscopic sur- 
veys in the mid-1980' s determined warm ISM detection frequencies in the cores of ear ly-type 
galaxies of 30%-50% jOaldwelllll98~l Iphillips et~al~1ll986l : Iveron-Cettv fc Veronlll986h . Sub- 
sequent narrow-band Ha imaging studies revealed extended emission in general agreement 
with the the stellar morphology and d etection rates consistent with the 50% ra te determined 
from the early spectroscopic surveys (jDemoulin-Ulrich et al.lll984j ; lKimlll989l ). 



1.2. Key Warm ISM Studies 



There are three large sample, multi- wavelength studies of early-t ype galaxies that have 

provid ed a foundation for our understanding of the cold and warm ISM. iTrinchieri fc di Serego Alighieri 
fll99jj , hereafter TdSA) conducted one of the first quantitative early-type galaxy emission- 
line studies with Ha imaging an d long-slit optical spectr osco py of 13 early-type galaxies 
known to have hot X-ray halos. ( Goudfrooij et al. 1994a b c]; Goudfrooij &: de Jong 1995 . 
hereafter G94-5) investigated the origin and evolution of the ISM in an extensive study of a 
complete sample of 56 elliptical galaxies using broad-band imaging, narrow-band Ha imag- 
ing, long slit spectroscopy, and Infrared Astronomical Satellite (IRAS) far-IR imaging. A 
study of the warm and cold phases of the ISM in early-type galaxies was carried as part 



of th e ESQ Key-Program 1-004-43K (IMacchetto et al.lll996l : iFerrari et all 1 19991 : ICaon et al. 



20001 ; IFerrari et al.ll2002l . hereafter MFCF). Similar to G94-5, the MFCF study probed the 
origin and evolution of the warm and cold ISM phases through a combination of broad-band, 
narrow-band and far-IR imaging and long-slit spectroscopy of 73 early-type galaxies. The 
two later studies had the benefit of FIR imaging capabilities and added the cold ISM (i.e. 
dust) to the discussion of the ISM in early- type galaxies. These studies sought to determine 
the frequency of detection of warm and cold ISM in early-type galaxies and then measure 
fundamental properties, such as mass, density, and luminosity. By correlating these proper- 
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ties with stellar and hot ISM properties, inferences are made about the origin and evolution 
of the ISM and how it relates to the origin and evolution of the galaxy itself. Below we briefly 
summarize the findings of these key studies and note the areas where the studies agree and 
where questions remain. 

The three key studies find similar fractions of their samples that host a warm ISM, 
ranging from 60%-80%, despite the differences in sample selection criteria and sizes of the 
samples. TdSA found 75% of the 13 early-type galaxies surveyed had detectable Ha emission 
which was extended over a range of 5-10 kpc. The extent and elliptical shape of the emission 
confirmed the galaxy- wide nature of this phenomena and excluded a nuclear confined process, 
such as AGN activity. The TdSA survey was conducted on galaxies known to have significant 
masses of hot ISM, and the authors concede that the high detection rate is likely connected 
to sample selection. The G94-5 optically complete sample (chosen without regard for X- 
ray luminosity) offers a less biased measure of the frequency of the warm ISM; 60% of the 
early-type galaxies surveyed were detected to have a warm ISM, while a cold (T~ 10K) 
dust component is detected in 40% of these galaxies. Through the FIR imaging, the G94-5 
study was able to show that unlike spiral galaxies, early-type galaxies display a smooth dust 
component distributed throughout the galaxy. The MFCF authors report an ionized gas 
detection frequency at 78% for the 73 galaxies in the optically complete sample. The cold 
dust component was detected at a 75% detection rate with smooth, extended morphologies 
similar to that of the warm ISM. Two key observations link the warm and the cold ISM 
together; first, the similar spatial extent and morphology of the cold ISM (dust) and warm 
ISM and second, the fact that the presence of dust is a strong predictor of the presence of 
emission lines. It is concluded that these ISM phases are closely related and likely of the 
same origin. 

The mass of the warm ISM is orders of magnitude below the X-ray emitting hot ISM and 
the stellar mass. TdSA calculated a warm ISM gas mass by assuming case-B recombination, 
a filling factor of 10~ 5 — 10~ 6 , and densities corresponding to the pressure implied by the 
X-ray observations. The resultant masses were calculated to be ~ 10 5 — 1O 6 M for a typical 
galaxy in their sample. The Ha imaging in G94-5 revealed warm ISM masses similar to the 
FIR determined dust masses of ~ 10 4 — 10 5 M Q . The MFCF study determined masses for the 
warm ISM in the range of 10 3 — 10 5 M Q . Masses of dust at 10 5 M Q have direct implications 
on the origin of the dust as the dust destruction rate is high in the presence of the hot ISM 
and is further discussed below. 

The correlations between the luminosities of L FIR , L Ha , L x , L B and L X /L B explore the 
links between the difference phases of the ISM and the stars. For numerous of reasons, the 
luminosity correlations have been difficult to accurately measure and there are contradictions 
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between the key studies. TdSA observed that galaxies with more massive X-ray halos had 
stronger line emission, indicating that the warm gas is somehow linked with the hot gas. 
No other luminosity correlations were observed. In a related work to G9 4-5, the authors 



find a strong anti-correlation between the masses of dust and X-ray gas (jGoudfrooijl 11994 



1 9 9 5 ) . However in a later study with a larger sample (based in part on MFCF data), 
Goudfrooijl (119971 ) reveals a positive correlation between the warm ISM with both the hot 
ISM and the stellar blue luminosity. MFCF determined the luminosity of the ionized gas to 
be correlated with both blue and X-ray luminosity. The scatter in the Ln a and Lb relation 
is greatly reduced when the blue luminosity is confined to the same region as the emission 
lines. The scatter in all the to Lx / Lb relations is large and plagued by upper limits on 
non-detections. 

In general, the three studies agree that the ionization source for the warm ISM is 
consistent with PAGB stars, but all comment that the data are unable to exclude other 
excitation sources such as internal shocks, AGN heating or electron conduction from the hot 
gas. TdSA show that the observed lin e intens i ties ar e a factor of ~ 10 2 — 10 3 too large to result 



from cooling flow luminosities (e.g. iFabianl (1l99ll )) and conclude that UV flux from post 



asymptotic giant branch (PAGB) stars are the most likely ionizing mechanism. The MFCF 
observed correlation between Ha and blue luminosity is an indication of photoionization of 
the warm ISM from PAGB stars. However, the observed correlation with X-ray luminosity is 



consi stent with warm ISM excitation via electron conduction from the hot ISM (ISparks et al. 



19891 ); Neither ionization mechanism is preferred in MFCF's large sample of 73 galaxies. See 
Section f4. 1.11 for a continued discussion on the warm ISM ionization source. 

The key studies narrow down the origin of this dust and warm gas to either be mass 
shed from AGB stars or the result of small galaxy accretion and mergers. G94-5 reasons 
that cold dust originated in the circumstellar envelopes of AGB stars are quickly sputtered 
away by the hot IS M with a typical dust gr ain lifetime (~ 10 7 yrs) much shorter than 
the age of the galaxy (jDraine fc Salpeterlll979l ). And further, the balance between the dust 
i njection and destruction r ate does not lead to the dust masses inferred by the far-IR imaging 
( IFaber &: Gallagher! Il976l ; iKnapp et al.l Il992l ) . The surprising result from the confluence of 
the G94-5 studies is that the majority of the dust in elliptical galaxies is external in origin. 
Since the distribution of the ionized gas (i.e. the warm ISM) is observed to be similar to the 
dust morphology, it is presumed that the warm ISM is likewise external to the galaxy, and 
will only be seen in galaxies that either have recent merger activity or significant inflows. The 
high detection frequency of the warm ISM, implies mergers and smaller galaxy acquisitions 
are common in early-type galaxies. Adding to the puzzle, in a subsample of the gas-rich 
galaxies, the warm ISM was frequently observed to be kinematically disti nct from the mai n 
stellar population, also strongly suggesting an external origin for this gas (jCaon et al.ll2000l ). 
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1.3. Warm ISM Open Issues 

The warm ISM is important to study in early-type galaxies and despite the ambitious 
efforts of the studies described above and subsequent works, several fundamental questions 
remain. It is critical to resolve these open issues since the intertwined relationships between 
the cold, warm, and hot gas phases and the stellar population ultimately defines galaxy 
evolution. Even though the warm ISM mass is orders of magnitude less than the stellar 
mass, discriminating between internal versus external origin of this material provides pivotal 
data in formation models. Also, even though minute in mass compared to the cold and hot 
phases of the ISM, the optical line emissions of the warm ISM can locally dominate the 
energy output for a region. Finally, the optical emission lines of the warm ISM provides 
information on the heating and cooling in these systems which demands further exploration 
since the observed hot ISM properties along with the non-detection of cooling flows creates 
a "cooling crisis." 

The metallicity of the warm IS M is unknown . Beca use the metallicities of the hot 
ISM and the stellar metallicities (cf. Trager et al. ( 2000bl la])) are known, determining the 



metallicity of the warm ISM can create a link between the gas phases and the stellar popula- 
tion. Also, the metallicity of the warm ISM contains important discriminating information 
concerning the outstanding issues that the key warm ISM studies were not able to resolve. 
Specifically, the line ratios of the ionization species provides information on the excitation 
mechanism of the gas; currently, both electron conduction and photoionization remain viable 
excitation mechanism. 

Also the metallicity of the warm ISM contains clues to the internal versus external origin 
debate. The gas injected into the ISM from s tellar mass loss will g ener ally be mo r e met al 



rich than gas accreted from dwarf galaxies (cf. iLequeux et al.l ( 119791 ) and lAaronsonl (119861 )). 



Because of the important information contained in the metallicity of the warm ISM 
and the unresolved issues concerning its nature and origin, we have conducted a program 
to determine the oxygen metallicity of the T ~ 10 4 K gas in a small sample of early-type 
galaxies. Section [2] describes the sample and observations. The data processing techniques 
and spectral analysis is discussed in Section [31 Results from the emission line analysis is 
presented in Section HI In Section O we conclude with a discussion of the implications of 
this work on both the origin of the warm ISM and the relationships between the stellar 
population and the other phases of the ISM. 
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2. Observations 

The MDM 2.4 meter Hiltner Telescope coupled with a Boiler and Chivens spectrograph 
was used to observe 21 nearby, early- type galaxies in four observation runs over a two year 
period from 2000 to 20 02. The sample was select ed to span a range of hot ISM properties 



as determined from the iBrown fc Bregmanl (119981 ) X-ray study as well as a range of warm 



ISM properties determined in the narrow-band Ha imaging surveys by TdSA, G94-5 and 
MFCF. The galaxies with their relevant fundamental properties are listed in Table [TJ Our 
warm ISM isolation and detection method depends on observing galaxies that contain no 
emission line gas, therefore, several galaxies were chosen to have very weak detections of the 
hot and warm phases. 

The configuration of the spectrograph was chosen to maximize galaxy light input with 
a 2.1" wide by 5' long slit, while retaining the ability to discriminate between the Ha and 
[N II] A6583 emission lines. The available and appropriate grating for these requirements 
was a 350 lines/mm grating blazed at 4026 A, resulting in 7.1 A/pixel and a spectral range 
of ~ 1600A over the 1200x800 pixel Loral CCD. The CCD was characterized by relatively 
low read noise at seven electrons with the nominal gain set at 2.1 electrons per ADU. The 
one drawback to the chosen spectrograph configuration is that it was necessary to obtain 
separate blue and red spectra for each galaxy in order to obtain all of the important emission 
lines from [O II]A3727 to [Si II]AA6717, 6731. The grating is servo controlled from the con- 
trol computer and our tests indicated a 5A accuracy in repositioning. Therefore, internal 
HgNe and Ne lamps were observed before and after each grating reposition for wavelength 
calibration. Additional calibrations included bias frames, evening and morning twilight flats 
when the skies were clear, internal flats illuminated from an incandescent source, and spec- 
trophotometric standards. 

The four observation runs were in March 2000, May 2000, February 2001, and February 
2002. Galaxy exposure times and program-type (E = emission-line, T = Template, or L 
= LINER) are listed in Table [2J The meaning of the program types is discussed below. 
During the first three runs, the nights were cloudy and over half of the full run was lost 
to weather and instrument problems. Some of these observations are of marginal value and 
only reported here for completeness. The final run was entirely photometric. 

The observing strategy was chosen to suppress spurious emission line detections and ob- 
tain uniform observations of program and template galaxies. For each galaxy, we attempted 
to obtain five integrations, ensuring that reasonable statistics could be employed to eliminate 
cosmic rays from the combined data. Because the emission lines are faint and spread over 
only a few pixels in both spectral and spatial directions, a fortuitous single cosmic ray can 
result in a spurious warm ISM detection. The total integration time was chosen to result in a 
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Table 1. Early- type Galaxy Properties: Warm Phase Sample 



Galaxy 




RA a 


Dec 3 






D b 


logL B /L Q h 


log L x c 


logL Ha d 


H a Reference c 




(J2000.0) 


(J2000.0) 




(km/s) 


log (erg s x ) 






NGC 1407 


03 


40 


11.8 


-18 34 48 


10.57 


1990 ± 187 


11.16 


41.34 


39.32 


M96 


NGC 1600 


04 


31 


39.8 


-05 05 10 


11.79 


4019 ± 489 


10.67 


40.84 1 


40.15 


M96 


NGC 2768 


09 


11 


37.5 


+60 02 15 


10.93 


1532 ± 325 


10.79 


40.41 






NGC 3115 


10 


05 


13.9 


-07 43 07 


9.95 


1021 + 215 


10.83 


39.74 


< 37.6 


M96 


NGC 3377 


10 


47 


42.3 


+13 


59 


08 


11.13 


857+ 126 


10.21 


39.42 


38.95 


G94 


NGC 3379 


10 


17 


49.6 


+12 


34 


55 


10.43 


857+ 126 


10.49 


39.78 


39.04 


M96 


NGC 3489 


11 


00 


18.3 


+13 


54 


05 


11. 12 2 


1039 ± 101 2 


10.35 2 




39.34 


M96 


NGC 3607 


11 


16 


54.3 


+18 


03 


10 


10.53 


1991 ± 242 


11.18 


40.82 


39.92 


M96 


NGC 4125 


12 


08 


05.8 


+65 


10 


27 


10.58 


1986 ± 295 


11.16 


41.01 


40.30 


T91 


NGC 4261 


12 


42 


02.4 


+11 


38 


48 


11.32 


2783 ± 590 


10.35 


41. 18 2 


39.38 


G94 


NGC 4374 


12 


25 


03.7 


+12 


53 


13 


10.13 


1333 ± 71 


10.99 


41.09 


39.56 


G94 


NGC 4406 


12 


26 


11.7 


+12 


56 


16 


9.87 


1333 ± 71 


11.10 


41.80 


40.50 


T91 


NGC 4472 


12 


29 


46.8 


+08 


00 


02 


9.32 


1333 ± 71 


11.32 


41.77 


39.60 


T91 


NGC 4494 


12 


31 


24.1 


+25 


46 


28 


10.69 


695 ± 147 


10.20 


39.28 






NGC 4552 


12 


35 


39.8 


+12 


33 


23 


10.84 


1333 ± 71 


10.71 


40.92 


39.26 


M96 


NGC 4636 


12 


12 


50.0 


+02 


41 


17 


10.20 


1333 ± 71 


10.96 


41.81 


39.69 


M96 


NGC 4649 


12 


43 


39.6 


+11 


33 


09 


9.77 


1333 ± 71 


10.96 


41.48 


39.83 


T91 


NGC 4697 


12 


18 


35.9 


-05 48 02 


10.03 


794 ± 168 


10.58 


40.13 


39.63 


G94 


NGC 5044 


13 


15 


23.9 


-16 23 08 


11.25 


2982 ± 314 


10.34 


42.39 


40.73 


M96 


NGC 5322 


13 


49 


15.2 


+60 11 26 


11.09 


1661 ± 352 


10.80 


40.11 


39.74 


G94 


NGC 5846 


15 


06 


29.2 


+01 36 21 


10.67 


2336 ± 284 


11.26 


42.01 


40.25 


M96 



a Values taken from NED (NASA/IPAC Extragalactic Database). 

b Radial velocity distance corrected for local flows from lFaber et al . (1989) unless otherwise noted. 
c X-ray luminosity from lBrown fc Bregmar] l|l998h unless otherwise noted. 
d Hct narrow-band imaging luminosity. 

e Re ference for Hq narrow-band ima ging luminos ity with the following abbrevia t ions: M96 = Macchctto et all 
l ll996h . G94 = lGoudfrooii et ail l ll994bl '), and T91 = lTfinchieri fc di Serego Alighieril Jl99lh . 



lO'Sullivan et al.l ll200ll) . 
2 RC3. lde Vaucouleurs et all dl99ll l. 



- 8- 



Table 2. MDM 2.4 Meter Observations 



Galaxy 


Exp Blue (s) 


Exp Red (s) 


Emission Type a 


NGC 1407 


8400 


6600 


T 


NGC 1600 


8400 


5400 


W 


NGC 2768 


12000 


11400 


L 


NGC 3115 


6000 


4500 


T 


NGC 3377 


7200 


6000 


W 


NGC 3379 


6000 


6000 


T 


NGC 3489 


6000 


4500 


E 


NGC 3607 


9600 


7200 


E 


NGC 4125 


12000 


10800 


L 


NGC 4261 


6000 


4500 


E 


NGC 4374 


6000 


4500 


E 


NGC 4406 


9600 


6000 


W 


NGC 4472 


9600 


6000 


W 


NGC 4494 


1200 





w 


NGC 4552 


6000 


6600 


w 


NGC 4636 


6000 


9900 


E 


NGC 4649 


6000 


4500 


W 


NGC 4697 


6000 


1800 


W 


NGC 5044 


6000 


3600 


E 


NGC 5322 


9600 


9900 


W 


NGC 5846 


7200 


6600 


E 



a Emission Type. Column display how the observations 
were ultimately used: T= template galaxy, L=LINER galaxy, 
E=cmission galaxy, W=weak emission. 
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S/N of 5-10 in the template subtracted spectra (See Section [3]). We did not constrain the an- 
gle of the slit due to complications with the instrument rotator during the first two runs; this 
has little programatic impact because of the choice of a wide slit, and the non-preferential 
orientation of the warm ISM observed in the narrow-band imaging surveys (TdSA, G94-5, 
MFCF). 



3. Reductions and Spectral Analysis 

The data were reduced in the standard manner using tasks within IRAF (overscan 
fitting and subtraction, bias frame construction and subtraction). Each of the runs was 
calibrated separately but in a similar manner. Internal lamps were used to create a response 
image and divided through the data, eliminating differences in pixel-to-pixel sensitivity on 
the CCD. Twilight flats were used to construct an illumination image, correcting for large 
scale structure and slit illumination. For the first two runs, the incandescent internal lamp 
produced too few photons in the far blue and only added noise to the data. For these runs 
we used wavelength calibrated twilight flats to produce response images. A HgNe lamp was 
used for wavelength calibration of blue-tuned observations and a Ne lamp was used for red- 
tuned observations. Residuals in the linear, spectral solution were much better (typically 
< 1A) than the resolution element (~ 8A) over the entire chip. 

The galaxies were summed over the spatial dimension based on the visual inspection 
of the 2-D spectra, with background regions selected from the outer regions of the chip. 
Typically the summed regions for the galaxies were about 1.2' of the 5.2' unvignetted slit 
l ength. The spectra were flux calibrated from observations of spectrophotometric standards 



(IMassey et al.lll988l ) and corrected for extinction with the standard correction for Kitt Peak 



distributed with IRAF 2.11. 

Custom code was developed to implement the combining of the individual 1-D spectra, 
template matching and subtraction, as well as the extraction of line fluxes. The 1-D spectra 
combination procedure employed an iterative 1.5a rejection algorithm which eliminated cos- 
mic rays that were overlooked in the aperture summing rejection methods. The wavelength 
calibration was tuned to night sky lines and the measured redshift was checked and tuned to 
strong stellar population lines (e.g. Ca H&K in the blue and [N II] A6583 in galaxies showing 
emission lines in the red). An example summed spectra is shown in Figured! 

Although we have upper-limits from the Ha imaging surveys, it is unknown which 
galaxies are void of warm ISM emission. Therefore, the determination of template galaxies, 
or galaxies without any detected warm ISM, is an involved and iterative process. First the 
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[N II] 




4000 5000 6000 7000 

Wavelength (A) 

Fig. 1. — Composite Blue and Red Spectra of NGC3489. Important emission lines that are 
visible in the un-subtracted data are marked aobve the continuum level. In addition, calcium 
H & K lines are marked for reference. 
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Ha and [N II] A6583 region of the spectra is examined for obvious emission. Using the local 
continuum, upper limits on the line fluxes are determined for galaxies with no obvious Ha 
and [N II] A6583 emission. Ten of the original 21 galaxies match this criteria and act as 
potential template galaxies. We take these ten galaxies and combine them into a master 
template galaxy. When combined into the template galaxy, the strongest of these weak 
emitters will not be included in the template creation because of the aggressive clipping 
algorithm in the combination routine. To determine the the weak emitters the template is 
subtracted from the individual galaxy observations and examined for emission lines. Three of 
the ten template galaxies revealed oxygen lines and were removed from the master template. 
This process was repeated with the emitters removed and reveal four more weak emitters, 
leaving three galaxies without emission line activity: NGC 1407, NGC 3115, NGC 3379. 
As an independent check of our template definition methodology, all seven of the potential 
template galaxies were subtracted from the galaxy with the strongest emission lines (NGC 
2768) and we determined that the three template galaxies selected resulted in the high 
reported fluxes (i.e. lowest background). 

The three template galaxies, NGC 1407, NGC 3115 and NGC 3379, lack of warm emis- 
sion is consistent with results from the narrow-band imaging detections for these galaxies. 
NGC 3379 is reported to have a weak Ha in MFCF but an upper limit below this reported 
detection is determined in G94-5. NGC 3115 is detected only as an upper limit by MFCF 
and has been previously used as a stellar t emplate galaxy by the Palomar spectral survey of 



the nuclear regions of 500 nearby galaxies (IHo et al.lll993l ). NGC 1407 is marginally detected 



in both the G94-5 and MFCF studies, but the fluxes reported differ by a factor of five. The 
inconsistencies of weak detections in the narrow-band surveys reveal intrinsic limitations of 
the narrow-band technique. Broad-band red imaging is used to subtract an underlying stel- 
lar population from the narrow-band data; a process which is susceptible to scaling errors. 
The stellar template constructed from our sample contains < 10~ 16 erg s^ 1 cmT 2 combined 
Ha and [N II] A6583 flux; this is an order of magnitude more sensitive as a non-detection 
than the narrow-band imaging surveys. 

Once the template galaxies are defined, the stellar population was subtracted from each 
of the other galaxies, revealing emission lines and any differences in stellar populations. An 
example subtraction for NGC 4374 is shown in Figure [2J In all galaxies, the residuals in 
the subtractions are flat with large-scale variation at a 8-10% level in the blue and better 
than 5% in the red. Once the subtractions are made, the relative line fluxes are measured 
by defining a low order polynomial to describe the local background and using a gaussian 
profile to extract the line flux. A gaussian decomposition is only necessary in the crowded [N 
II]A6548, Ha, [N II] A6583 region and the fitting method gives equivalent results to raw counts 
above a background for all other lines. Note that because of the coarse spectral sampling and 
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the relatively low velocity dispersions, the width of the gaussian profiles contain no useful 
kinematic information and is thus not reported. In eight of the non-template galaxies, the 
signal-to-noise in the data are insufficient for emission line work, or only Ha is observed. 
For the rest of the sample, extracted relative line fluxes normalized to H/3 are reported 
in Table [3j The repo rted flux ratios have bee n de reddened based on a n average value 
for the reddening from iBurstein fc Heilesl (119841 ) and ISchlegel et al.l (119981 ). The observed 
relative fluxes reported here do not correct for internal galaxy reddening, displaying the 
Balmer decrements as observed; Below we apply the correction (See Section |4~21) . The errors 
reported are determined by calculating the range of acceptable background levels within the 
noise of the local continuum for isolated lines and combining this background uncertainty in 
quadrature with the range of acceptable gaussian widths that reproduce the total observed 
line flux for the crowded [N \V\-Ha complex. 



Warm ISM Emission Line Results 



4.1. Excitation Mechanism 



4-1.1. Normal Galaxies 



There are many possible ionization mechanisms for the warm ISM that have been pre- 
viously explored and ruled out by th e key warm ISM studies (Section 11.21) and other works. 
Photoionization by hot young stars (jKiml Il989l ; IShieldd Il99ll ) is inconsistent with both the 



optic al colors (G94-5; MFCF) and integrated optical sp ectra ( Sadler! 119871; Heckman et al. 



19881 1 . Ionization from radiation associated with an AGN (IFosbury et al.lll982l ) is ruled out by 
the large radial extent of the warm ISM observed in the narrow-band imaging (TdSA; G94- 



Table 3. Observed Nebular Emission Line Ratios Relative to Hj3 



Galaxy 


[O II] 
A3727 


[Ne III] 

A3869 


Hfi 
A4862 


[O III] 

A4959 


[O III] 
A5007 


[O I] 
A6300 


[N II] 
A6548 


Ha 
A6563 


[N II] 

A6583 


[S II] 
AA6717, 6731 


NGC 2768 


70 


± 


7 


14 ±6 


10 ± 


3 


1±4 


24 ±4 


8±2 


20 ±4 


29 ±5 


45 


± 


7 


50 ±8 


NGC 3489 


19 


± 


5 


2 ± 3 


10 ± 


2 


2 ± 5 


13 ±5 


3±2 


30 ±8 


34 ±5 


57 


± 


8 


37 ±7 


NGC 3607 


32 


± 


8 


6 ± 3 


10 ± 


3 


7±3 


20 ±4 


2±4 


25 ±5 


48 ±6 


59 


± 


8 


22 ± 5 


NGC 4125 


74 


± 


8 


14 ±4 


10 ± 


2 


23 ± 5 


71 ± 10 


3±3 


19 ±4 


34 ±8 


40 


± 


5 


23 ± 7 


NGC 4261 


16 


± 


4 


6 ± 5 


10 ± 


3 


6 ± 2 


12 ±3 


3±2 


15 ±5 


22 ±8 


39 


± 


5 


15 ± 6 


NGC 4374 


16 


± 


4 


< 1 


10 ± 


2 


5 ± 2 


12 ±4 


2±4 


21 ±9 


35 ±6 


49 


± 


8 


29 ± 5 


NGC 4636 


24 


± 


5 


2±4 


10 ± 


3 


4±4 


13 ±3 


2±5 


10 ±3 


33 ±5 


49 


± 


7 


21 ± 5 


NGC 5044 


24 


± 


3 


2 ± 3 


10 ± 


2 


6 ± 2 


20 ±3 


3±2 


33 ±8 


32 ±6 


53 


± 


8 


33 ± 5 


NGC 5846 


26 


± 


5 


2 ± 3 


10 ± 


4 


6 ± 3 


14 ±4 


< 1 


8±5 


29 ±4 


36 


± 


8 


13 ± 10 
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5; MFCF). Shock excitation (IFord fc Butcherl Il979l ; iHeckman et al.l Il989l ) would produce 
filamentary ionization regions, where smooth emission is observed (TdSA; G94-5; MFCF). 
Further, the the input energy from sho cks is too low by tw o orders of magnitude for typ- 
ical densities and velocities as noted by ISparks et al.l (119891 ) . Condensation out of the hot 
phase and into warm phase filame nts has been investigated for a number of X-ray clusters 



flCowie et alJll980 



Hu et al 



1985) but th e large number of recombinations required per hy- 



Johnstone et al.l (119871 ) ; TdSA and refe rences therein) is in consistent with 



drogen atom (e.g. 

the observed recombination rates for early-type galaxies (ISparks et al.lll989l ). 



The three remaining viable excitation mechanisms for the war m ISM are ph o toion ization 



from PAGB stars, electron conduction from the hot phase (e.g. ISparks et al.l (1989V), and 



photoionization by ext reme UV (EUV) photons from the hot phase ( IVoit fc Donahue! 11990 



Donahue &: Voitl Il99ll ) . The data in the key warm ISM studies are equally consistent with 
all of these ionization mechanisms. However, photoionization by PAGB stars is the preferred 
mechanism since these old stars are known to populate the galaxies and given the signifi- 
cant correlation of Lu a with Lb within a given region (MFCF) strongly indicates a stellar 
ionizing source. Further, the presence of the far ultraviolet (UV) flux has been confirmed 
i n a small sample of 30 early-type galaxies with International Ultraviolet Expl orer (IUE) 
(Bur stein et al.l (119881 ); For a review of Far UV flux in early-type galaxies see lO'Connell 



(|1999|)). 



However, to adopt PAGB stars as the ioni zing source, exclu des the observed correlations 



between Ha and X-ray luminosity (TdSA; (IGoudfrooiil 11997), MFCF). The correlations 



between the X-ray and warm ISM morphologies (jSingh et al. 



1995; Trinchieri et al. 1997 



Trinchieri fc Goudfrooijl |2002| ) provide additional support for photoionization of the warm 
ISM by the EUV flux from the hot ISM, but it is possible to construct physical conditions 
that reproduce these morphological links (e.g. cooling of the hot ISM). Although there 
have been detected correlations between the L Ha and L X /L B , the correlations contain much 
more scatter and are less significant than the correlation between L# Q with Lb within a given 
region. The large amount of scatter in the X-ray relationship and the known high Lx/Lb 
galaxies which contain weak or no Ha emission and strong Ha emitters with low Lx/Lb 
provide insight that the hot and warm ISM link is more tenuous than the link between 
the st ars and the war m ISM. Finally, the key studies in in Section 11.21 and subsequent 
work (jSarzi et al.ll2006l ). indicate a strong connection between the dust (cold ISM) and the 
emission line regions, but there are not strong correl ations between L^ jr and Lx/Lb in any 
of the key studies or in recent Spitzer Observations (ITemi et al.l 120071 ) . 



Emission line ratios of the different ionization species have the ability to discriminate 
between collisionally excitation and photoionization. Comparing our sample's line ratios to 
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line ratios predict e d for photoionization models of I Allen et al.l (19981) . lAUen et al.l (119991 ). 
and iBinette et al.l (119961 ) and shock ionization of lAUen et al.l (120081 ). we conclude the data 
are consistent with photoionization. Unfortunately, these data are not discriminating due to 
the large errors in the data and the region populated by these galaxies happens to land in 
an overlap region for collisional and photoionization. 

Additional support f or the photoionizatio n of this gas is seen with the Far Ultraviolet 
Space Explorer (FUSE) (IBregman et al.ll2005l ). The observed UV radiation field produces 
sufficient photons shortward of 9 12 A to ionize the warm ISM and is consistent with the 
number of predicted extreme horizontal-branch stars, or equivalently PAGB stars. 

When the data presented here is combined with the literature, it presents a reasonable 
case for photoionization of the warm ISM. This is the preferred ionization mechanism in the 
literature and in discriminating between the two sources of photoionization, PAGB stars or 
the hot ISM, it would seem more probable that PAGB stars provide the necessary ionizing 
photons. 



4.1.2. LINER Galaxies 



Four of the galaxies listed in Table [3] are reported as low ionization nucl ear emission- 
l ine regions (LI NERs) galaxies in the literature: NGC 2768 ([lieckman! [l98oh . NGC 3115 
(THo et al.l[l993h . and NGC 4125 and NGC 4261 JHo et alJll995h . The question of identical 
treatment of these LINER galaxies and the normal galaxies arises. LINERs are strictly 
defined to be galaxies with line ratios of [ O II] A3727 great er than [O III] A5007 and [O 
I]A6100 greater than one-third [O III] A5007 (IHeckmanl Il980l ) . While useful for separating 
classes of objects, the LINER definition as noted by its creator, Heckman, has no physical 
basis and in reality only provides a crude boundary between photoionization and a low-level 
AGN activity. 

For our purposes, the LINER identification acts as a flag for further investigation, and 
not as an automatic exclusion filter. The line ratios for the four LINER galaxies are consistent 
within the photoionization, but the excitation mechanism and thus the valid regions in 
diagnostic diagrams are still uncertain a nd being debated for L INERs (for a review of the 
physical mechanisms behind LINERs see iFilippenkol (119961 ) and iBarthl (120021 ) ) . Although it 
is not entirely clear what to expect for LINER galaxies, the presence of strong collisionally 
excited lines casts a cloud of doubt of pure photoionization. The [Ne III] A3869 line is an 
indicator of the importance of collisional excitation because the collisional cross-section is 
over an order of magnitude larger than the photoionization transition probability for this 
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line. NGC 2768 and NGC 4125 have significant [Ne III]A3869 fluxes and thus we conclude 
that these systems have a more complex physical structure and excitation mechanism than 
photoionization. These two galaxies are excluded in the subsequent sections where pure 
photoionization is assumed. Note that NGC 4125 and NGC 4261 show no distinguishing 
line ratios from the normal galaxies and, the data presented here do not support the LINER 
classification. These two galaxies were classified as LINERs by ( Ho et al.lll995l ) which had 
a similar observational aperture as this study, but only had wavelength coverage down to 
4200A, excluding the useful [O II]A3727 information. 



4.2. Oxygen Metallicity 

For those galaxies where photoionization is the dominant ionization mechanism, the 
line ratios can be used to determine the metallicity of the warm phase gas. There are two 
critical assumptions to be addressed before using these line ratios and photoionization models 
to determine global metallicities. The first assumption is that UV flux from PAGB stars 
ionizing a diffuse ISM on a galactic scale is well represented by photoionization models that 
were developed for Stromgren-sp here- type H II regions. First, we note that the models are 
determining the balance between photoionization by UV photons and recombinations to the 
ions; this is fundamentally the same equation in a galaxy-wide ISM as an HII region, given a 
source of ionizing photons internal to a volume of gas. The UV flux and electron densities of 
the two environments are similar, a lthough, PAGB s tars have different spectral shapes than 



high mass main sequence stars, but iMcGaughl (119911 ) shows that the photoionization models 



are largely insensitive to the spectral shape of the ionizing star. The models assume the gas 
is ejected from the ionizing star (and thus of the same metallicity), this is a feature used to 
calibrate the observed line ratios, but a direct link between the ionizing source and the ejecta 
is not a necessary condition to employ the models. We argue that the photoionization models 
will be applicable to the warm ISM, but we readily acknowledge that there are differences in 
these two environments which will lead to will lead to greater uncertainties in its application 
here. 

The second key assumption is that the analysis of the linear combination of many in- 
dividual diffuse ionized regions observed in an integrated galaxy spectra produces a result 



that is indicative of the mean properties of the observed regions. iKobulnicky et al.l (119991 ) 
conducted a study addressing these issues in nearby galaxies, observing H II regions in- 
dividually through small apertures and globally with drift scanning with a long slit. The 
authors conclude that global metallicity determinations, accurate to 0.2 dex, are possible 
with long-slit, integrated spectra. 
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Fig. 2. — Blue (left) and Red (right) spectra of NGC4374 with template galaxies subtracted, 
removing the stellar continuum. Large residuals at 5500A and 6800A are due different 
night sky lines positions in the rest wavelength of NGC4374 and template galaxies. The 
large residual at 3900A is due to the imperfect matching between the Ca H&K lines of the 
template galaxies and NGC 4374. 



Table 4. R23 parameter and Oxygen Metallicity 



Galaxy log R23 [O/H] Z 0x yge n /Z e 



NGC 3489 


0.63t° ;H 


-3.26 


± 


0.2 


1 10+0.49 
± - lo -0.30 


NGC 3607 


0.93iH? 


-3.66 


± 


0.2 


n 4 o+0.29 


NGC 4261 


51 + 06 
u - o± -0.07 


-3.14 


± 


0.2 


1 44+0-47 


NGC 4374 


0.58t° ;?g 


-3.20 


± 


0.2 


i-26l8:I! 


NGC 4636 


0.68t°;°? 


-3.31 


± 


0.2 


0.98t°;« 


NGC 5044 


7^+0.05 
u - ,o -0.06 


-3.37 


± 


0.2 


0.86l°i« 


NGC 5846 


n fis+ 08 

U - b °-0.09 


-3.30 


± 


0.2 


1 00+ ' 37 
l.UU_ 26 
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Metallicity determinations of ionized gas in individual H II regions are built upon pho- 
toionization models, determining ratios of ionized species relative to hydrogen. The most 
common d iagnostic is the R23 parameterizat ion, R23 = {{\Q I I) X3727+[0 III] AA4959, 5007)/H[3), 
defined by lPagel et all (119791 ) and refined by lMcGaughl (ll99ll ). This parameterization is con- 
venient since it employs the readily available optical oxygen lines, however, it requires the 
[O II]A3727/[N II] A6583 information to resolve a degeneracy present in the R23 models. 

We take the line ratios reported in Table p and correct for redden ing internal to the 
galaxy using the interstellar extinction cu rve of ISavage fc Mathid (Il979l ) and assuming case 
B ratio for the Ha to H(3 at T = 10 4 rtOsterbrockl Il989h . For NGC 4261 no reddening 
correction is applied, as the Balmer decrement indicates a non-physical negative reddening. 

These dereddened lines are used to generate R23 values for the galaxies. Because of 
the faint nature of the [O III] A4959 line, we used the conical, equilibrium value of 1/3 the 
line flux of [O III] A5007 when this line was buried in the noise of the subtracted spectra. 
All of the galaxies end up on the upper branch of the R23 curve as in dicated by t he [O 
II]A3727/[NII]A6583 diagnostic . To determine metalli cities, we employ iMcGaughl (1199ll ) 
calibration as parameterized by iKobulnicky et all (119991 ) and reported the metallicity results 
in Table HI For the zero point of t he solar metallicity we employ a sola r photospheric 
oxygen abundance of [0/H]= —3.3 (lAsplund et all 120051 ; IScott et all 120081 ). but note the 
uncertainty and current debate about the solar determinations. The errors reported on the 
oxygen metallicities in these galaxies are at least 0.2 dex as these are the dominating and 
intrinsic uncertainti es in the models u s ed to create the R23 grid and its application to global 
galaxy spectra (cf. IKobulnicky et all (119991 ) ) . The median metallicity for the seven early- 
type galaxies is O/O =1.O1. Within these 0.2 dex absolute errors, the none of the galaxies 
are significantly different from one another. 



5. Implications and Discussion 

The determined mean oxygen metallicity of solar for the warm phase of the ISM has 
implications for the origin of this gas and its relationship to other phases of the ISM and the 
stellar population. The limitations of this study are apparent with a small sample of seven 
galaxies, large 0.2 dex errors, and a less than ideal application of H II region models to the 
warm ISM. However, this work does present new information as it appears that there are 
no determination of the metallicity of the warm ISM in early-type galaxies in the literature. 



Sarzi et al.1 (120061 ) conducted the most recent large survey of ionized gas in early-type galaxies 
and conclude at the end of the study that the origin of the ionized gas is still unsolved and 
state that the metallicity provides an insightful clue. Below we consider the implications of 
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a near solar metallicity for the warm ISM. 



5.1. Origin of the Warm ISM 



A warm ISM with a solar metallicity argues against an external origin model for the 
emission line gas. The accretion galaxies involved in minor mergers are presumably dwarf 
irregular (dlrr) galaxies, since dwarf spheroidal (dSph) galaxies are most frequently devoid 
of gas. Gaseous and du sty galaxies with near solar metallicity are luminous (> 10 10 L o ) and 



massive ( IGarnettl 120021 ). Mergers of these types of objects produce profound disturbances, 
which are simply not observed in these early-type galaxies with Ha emission. It is possible 
that an intense gas enrichment by supernovae type-II occurs once low metallicity gas is 
accreted, however, the broad-band colors and global galaxy spectra (G94-5, MFCF) show 
little evidence for major star formation. Further, the observed su pernovae rates are too low 
to account for such a drastic enrichment (ICappellaro et al.lll999l ). 



We can compare a solar metallicity for the warm ISM with metallicities of the hot ISM 
and the stars. The stellar popula tion in early-type g alaxies is generally observed to have 
a solar or super-solar metallicity ( Trager et al. 2000bl Jah . A number of the seven galaxies 
observed in this work have stellar Fe and a metallicity measurements in the literature, 
although, there are some inconsistencies between authors (NGC 4261 is measured to have 
[FejH] of - 03, 0.275, 0.29, and 0.188 bvllrager et all (boOObh , iThomas et all (|2005h , Irlowell 
( 120051 ) . and lSanchez-Blazquez et al.1 (120061 ) . respectively.) Not all the galaxies have measured 
stellar metallicities but a five galaxy average from all of the determinations in the 2005 and 
2006 papers listed above is [Fe/H] = 0.22, ranging fro m \Fe/H] = . 099 to [Fe/H] = 0.44 
([a/H] is similar). It is important to remember that iTrager et al.1 (j2000al ) estimates the 
absolute error on stellar metallicity determinations close to a factor of two, while the relative 
errors are much smaller. The Chandra X-ray Observatory has the spatial resolution to resolve 
and remove the stellar component from the gas which is leading to bette r data and a revision 
of our (mis)understanding of the metallicity in the hot ISM. Recently, iHumphrey fe Buote 
( 120061 ) analyzed a sample of 28 early-type galaxies observed with Chandra and determine 
that the hot ISM abundances are consistent and perhaps higher than the stellar abundances. 
In summary, a 0.50 to 1.50 solar metallicity in the warm ISM is consistent with both hot 
ISM and stellar metallicity determinations. 

It has been thought that early-type galaxies would have a problem generating and 
retaining dust and warm gas, due to the hostile conditions provided by the hot ISM an d the 
low injection rate of new material into the ISM from stars. iParriott fe Bregmanl (120081 ) have 
shown that is it possible for dust to condense in the circumstellar envelopes of AGB stars 
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and then subsequently be injected into the hot ISM. Although this dust has a short lifetime 
in the hot ISM, the injection rate has been observed to be greater than the destruction rate 
and could account for the total observed mass given a simple evolution o f this injection rat e 



given that more massive stars evolved off the main sequence in the past (lAthey et al.ll2002l ). 
An internal origin for the dust is supported by the G94-5 and MFCF studies indicating 
that the majority (~ 90%) of the dust is distributed smoothly throughout the galaxy. The 
internal origin for the cold ISM is further confirmed by a 2.29/im CO absorption feature 
study that finds a lower veloci ty dispersion for the dus t than that of the stars in a sample of 



25 nearby early-type galaxies (jSilge fc Gebhardtll2003l ). Finally f rom the analysis o f Spitzer 



observations of 46 galaxies in the FIR (24/xm, 70/zm, and 160/xm). lTemi et al. J2007h propose 
an AGB mass loss, internal origin for the observed dust, which sinks to the core of the galaxy 
and is periodically driven out by AGN activity. This outward transport of dust may cool 
the hot gas down to the warm phase. 

The major stumbling block to a galaxy internal origin for dust and gas is the observations 
of gas kinematics which are distinct with respect to the stars. Recent observational evidence 
has suggested that there may be different evolutionary tracks for early-typ e galaxies. The 



classifi cation of early-type galaxies into fast and slow rotation categories by lEmsellem et al. 



( 120071 ) and the key ingredient of gas for the formation and evolution of fast rotators provides 
an important dividing line for future discussions on the origin of the ISM in early-type 
galaxies. 



The authors would like to acknowledge useful discussions with P. Goudfrooij and E. D. 
Miller. We would also like to thank an anonymous referee for suggestions that significantly 
improved this paper. 
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